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A NONLINEAR THEORY FOR AIRFOILS WITH 
TRAILING -EDGE JET FLAP 

Raymond L. Barger 
Langley Research Center 

SUMMARY 

A nonlinear procedure for computing the pressure distribution on an airfoil with a 
trailing-edge jet flap is described. The method is not restricted to thin airfoils or shal- 
low jet-deflection angles. Correlation with experiment indicates that the characteristics 
of the pressure distribution are predicted by the theory, but the effect of entrainment is 
overpredicted with the entrainment coefficient used. 

INTRODUCTION 

Interest in the jet flap stems from its effectiveness as a propulsive lift mechanism. 

A system, such as the jet flap, that increases the airfoil circulation is much more efficient 
than a pure jet reaction device. 

Theoretical studies of jet-flap airfoils have special importance because of the diffi- 
culties involved in wind-tunnel tests on such complex devices. The initial theory for an 
airfoil with a jet flap was developed by Spence (ref. 1) who used thin-airfoil calculations 
together with thickness corrections obtained from transformation theory. The effect of 
jet vorticity was calculated by distributing vortices along the X-axis (consistent with thin- 
airfoil theory). Comparison with experiment indicated better agreement for the lift 
coefficient than for the pressure distribution because of compensating errors in the upper 
and lower surface pressures. 

Spence did not treat the effects of jet entrainment, but such effects have been con- 
sidered by Wygnanski in reference 2. Wygnanski also used a thin-airfoil analysis, with 
the symmetric airfoil at zero incidence represented as a straight-line segment, and the 
jet by a line of sinks extending downstream from the airfoil on the X-axis. 

The present treatment includes the effects of both the jet deflection and the entrain- 
ment. The thin-airfoil assumptions are dropped so that arbitrary airfoil shapes and 
deflection angles are permitted. Boundary-layer effects are not included. 
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Fourier coefficients 

length parameter used in transformation of airfoil to approximate circle 

jet- momentum coefficient 

jet mass flow coefficient, — — 

pressure coefficient 

chord length (assumed unity in calculations) 
imaginary part of complex quantity 
imaginary unit 

jet slope at airfoil trailing edge, tan r 
index 

radius of circle into which airfoil is transformed 

real part of complex quantity 

radial coordinate in z'-plane 

sink strength in main region of jet 

sink strength associated with kth jet segment 

sink strength in mixing region of jet 

distance along jet 

local flow velocity in jet direction 

jet velocity 

complex velocity potential 
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conjugate of complex velocity 

value of w in circle plane excluding contribution due to supercirculation 
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total value of w in circle plane 

Cartesian coordinates in airfoil plane 

complex variable representing points in circle plane 

complex variable representing points in near-circle plane 

vortex strength representing supercirculation due to jet sheet 

strength of vortex located in kth segment of jet sheet 

vortex strength associated with camber and angle of attack 

vorticity 

jet- slot width 

complex variable representing points in airfoil plane 
angular coordinate in z 1 -plane 
curvature of jet sheet 
density 

empirical parameter used in computing entrainment coefficient 

jet-deflection angle 

angular coordinate in circle plane 

defined by relation r = ae 1 ^ 
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average value of ip 


Subscripts: 


airfoil plane 


circle plane 


due to entrainment 


free -stream conditions 


Bar over a symbol denotes complex conjugate. 


ANALYSIS 


Transformation of Airfoil Into Circle 


The first step in the analysis is to transform the airfoil into a circle by the method 
of Theodorsen (ref. 3). This method uses two consecutive mappings. The first is the 
Joukowski transformation 


which maps the airfoil in the £-plane into a contour approximating a circle in the z' -plane 
(fig. 1). The parameter a is determined so that if the airfoil is an ellipse of any thick- 
ness ratio, the mapped contour is an exact circle. (See ref. 3.) 

The complex variable z' is now written in polar coordinates 


and a new variable ip is defined by 


so that 


r = ae 1 ^ 


z’ = ae^ +i0 


n 


f 


Then the mapping of the contour to an exact circle in the z -plane (fig. 1) is accomplished 
by the relation 


z . = z eS( A n riB n)/ z " (2) 

where A n and B n are R n times the Fourier cosine and sine coefficients, respec- 
tively, of the function ip(8), which is obtained directly from the first (Joukowski) trans- 
formation. The radius R of the circle to which the airfoil is mapped is computed from 
the formula 


i 

1 


where if/ is the average value of i//. 

Reference 3 demonstrates that relations (1) and (2) map points on the airfoil into 
points on the circle in the z-plane. However, since equation (2) is not analytically 
invertible, points that are not on the airfoil are not easily mapped into the z-plane, and 
some iteration is required. The iteration procedure, however, is simple and rapidly 
convergent. 

On the other hand, the mapping of points from the circle plane to the airfoil plane 
is a straightforward process. Furthermore, velocities computed in the circle plane are 
simply converted to velocities in the airfoil plane by the relation 
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| dW ' 
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with 

dz 


given by differentiation of equations (1) and (2): 
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A n + iB n \ 
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Modeling the Jet Sheet 

To determine the shape of the jet and the circulation that results from the jet 
deflection, the jet is modeled as a vortex sheet. The derivation of the required formulas 
is given in reference 1. The known quantities are the jet -deflection angle t at the 
airfoil trailing edge and the jet-momentum coefficient 
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The shape of the vortex sheet and the distribution of vorticity along the sheet are 
to be determined so that 

(1) The boundary condition on the airfoil surface is satisfied. 

(2) The vortex sheet has the given deflection r at the trailing edge. 

(3) The curvature of the jet sheet is consistent with the distribution of vorticity 

U 2 C-/c 

along the sheet: that is, the relation y = — — 1— must be satisfied, where y and k 

2U 

are the local vorticity and curvature, respectively. 

An iterative procedure is required for calculating the jet shape and the vorticity 
distribution so as to meet all these conditions simultaneously. The jet sheet is approx- 
imated in the circle plane as a set of straight-line segments of equal length As c (see 
fig. 1). The first segment is determined so that the corresponding line in the airfoil 
plane has the given jet-deflection angle. For the first iteration, the remaining segments 
are chosen arbitrarily so that the slopes decrease exponentially with distance s along 
the jet. Then the vertex locations are determined, each point from the previous one, 
starting at the trailing edge. By moving the distance of one segment length in the direc- 
tion specified by the slope for that segment, the location of the next point is determined. 
The curvature distribution is also simply determined since the change in angle Ad at 
eacH vertex can be obtained from the slopes of the adjacent segments and the segment 

A f} 

lengths are equal. Thus, the curvature — can be calculated directly. A vortex is 
located at the center of each of the segments. 

Now consider the representation of the jet sheet in the airfoil plane. This jet 
sheet consists of an equal number of segments, but the segments are not of equal length. 
The first segment has the direction specified by the given value of r. A vortex is 
located on each segment with a strength determined by the approximation 


r a = r a A s a =o 


1U C 


2 A e„ 


2 Uo As, 


As a = i^: 
a 2 U a 


Ad, 


where A 0 a is the average change in angle at the vortex and is determined by the aver- 
age change in angle at the adjacent vertices. Since the mapping is conformal at each 
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vertex (although it is generally not conformal at the initial point located at the trailing 
edge), the change in angle is the same as the change in angle at the corresponding vertex 
in the circle plane. Thus, 


A0 a = A0 C 
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In general, it is not true that U a = U c , but for all cases calculated, U a differed 
so slightly from U c at the vortex locations that the adjustment in r because of the 
difference between U a and U c was essentially negligible. Therefore, the 
approximation 


r =l£« 

a 2 U r 


■A 9r 


(3) 


was adopted. This approximation shortens the calculation significantly as it permits the 
iteration to be carried out entirely in the circle plane. 

The iteration procedure requires the calculation of the velocity vector at each vor- 
tex location. The velocity equation in the circle plane is 
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(4) 


where the terms, in order, represent the undisturbed free stream, the doublet effect of 
the circle, the vortex due to the circulation associated with camber and angle of attack, 
the additional vortex associated with airfoil circulation resulting from the presence of 
the jet sheet, the vortices representing the jet sheet, and their images in the circle. 

Here, actually denotes the complex conjugate of the velocity variable. 

QZ 

On each iteration, the calculation requires several steps. First, the vortex 
strengths 1^ are determined by applying equation (3) at the kth segment. Then, 
since is known (see ref. 3), the only remaining parameter to be determined in the 
velocity equation is the "supercirculation" Ij. This quantity is obtained from the condi- 
tion that the first segment of the vortex sheet must have the direction specified by the 
given jet-deflection angle. The velocity at the center of the first segment can be written 


w T = w p + r jWq 


(5) 
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where w„ = — — , and w_ denotes all the remaining terms on the right side of equa- 

U x'7T l 7 * U 

tion (4). Since w T is the complex conjugate of the velocity variable, the relation 


Im(w T > 

m 1 = tan t = 

1 Re ( W TJ 


( 6 ) 


must be satisfied. Substituting from equation (5) into equation (6) and solving for r\ 
yields 


Inn w ' + m* Re-'w \ 

P 1 P' 

Im Wqj + m^Re Wp^ 


t 


With this value of Ij, the velocity vector can be computed at the center of each 
segment of the jet sheet; thereby the slope of each segment is fixed. This distribution 
of slopes is summed to obtain the jet shape for the next iteration, and it is differenced 
to obtain the curvature. Then, the curvature and the magnitude of the velocity are used 
to obtain the distribution of vortex strengths for the next iteration by means of 
equation (3). 

Considerable care is required in setting up the iterative calculation because of the 
dependence of the velocity on the curvature of the jet sheet. Any slight "kink" in the 
sheet arising in the calculation denotes an abnormally high local curvature and thus 
tends to trigger an instability leading to divergence. Approximating the sheet with poly- 
nomial or trigonometric series is not practical because increasing the accuracy of the 
approximation in the least-squares sense requires increasing the number of terms, 
thereby increasing the number of "wiggles" in the shape, and causing large local values 
of curvature. A spline fit generally results in even larger curvature variations. 

The smoothing process that was adopted used a smooth analytic shape on the first 
iteration. The results of each subsequent iteration were averaged with those of the 
previous iteration. This procedure proved satisfactory for values of Cj up to three. 
For larger values of Cj, it was necessary to start at three and gradually increase Cj 
during the iteration process. The convergence criterion adopted was that the value 
of Ij must vary less than 0.5 percent between successive iterations. The number of 
iterative cycles required has varied from 6 to 16 in the examples studied. 

The problem of convergence of the jet shape would probably be reduced if the vor- 
ticity were distributed over each segment rather than concentrated at a single vortex at 
the midpoint of the segment. 
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Velocity Due to Jet-Entrainment Effect 


After the jet shape has been established, one can compute the entrainment effect 
by distributing sinks along the jet. Formulas for the sink strengths per unit length, as 
given in reference 2, are 


S m = 0.064(V - Uco ) 

for the initial mixing region of the jet, and 


Sj = u„ 


3 (°j - c m) 


8cr(6/c) 


1/2 


( 7 ) 


( 8 ) 


for the main region of the jet, where the velocity profile is self-similar. These for- 
mulas are derived for a two-dimensional jet in a uniform free stream parallel to the jet 
axis, but some experimental results quoted in reference 2 indicate that they also repre- 
sent a good approximation for an inclined jet. The parameter ct is estimated to be 7.7 
from experimental data (ref. 2). Since formulas (7) and (8) give the sink strength per 
unit distance along the jet, the actual strengths are determined by multiplying these 
values by the corresponding distances along the jet. A sink is located on each of the 
straight-line segments, and the strengths are determined in the airfoil plane. For each 
of these sinks, there is a sink of equal strength at the corresponding point in the circle 
plane (ref. 4, sec. 8.50). In the circle plane, the velocity due to entrainment can be com- 
puted by the circle theorem which gives 



Thus, the boundary condition is satisfied by including, in addition to each sink, the image 
of the sink in the circle and a source of equal strength at the origin. The velocity due to 
entrainment is computed on the circle and added to the previously computed velocity 
distribution. 

It may be noted that the attempt to include the entrainment calculation in the 
iterative procedure for the jet shape is not practical because the jet sheet has to be 
represented as a distribution of vortices to obtain its shape, but it must be represented 
as a distribution of sinks to compute the entrainment effect. 
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Examples and Discussion 

The potential advantage to be gained by the use of a jet flap is illustrated by the 
example shown in figure 2. For this example, the thickness ratio is 0.13, Cj = 2.0, 
and r = 45°. One can see from the figure that the supercirculation due to the jet flap 
generates a large increase in lift. The effect of entrainment on the lower surface 
pressures tends to reduce the lift, but this effect is more than compensated for on the 
upper surface. The total increase in lift coefficient is about 500 percent, from 0.5 to 
approximately 2.5 percent. 

The only example of a measured pressure distribution given in reference 1 is for 
a 12.5-percent thick ellipse with a jet-deflection angle of 31.4°. When the nonlinear 
theory without entrainment is compared with Spence's modified thin-airfoil theory for 
this case (fig. 3), the results of the two theories compare very closely over most of the 
airfoil. This similarity was to be expected inasmuch as the jet-deflection angle was 
small, the airfoil was relatively thin, and thickness correction factors were applied to 
the thin-airfoil results. The nonlinear theory is somewhat closer to the experimental 
results in the regions of high curvature near the trailing edge. For a more conventional 
airfoil shape, the Spence theory would compare even more closely if accurate thickness 
corrections could be made. However, since the full transformation is required to obtain 
accurate thickness corrections for a general airfoil shape, the entire nonlinear theory 
can be used with little more trouble. The advantages of the nonlinear theory would also 
increase with increasing thickness ratio and jet-deflection angle. 

When the velocity due to entrainment is included in the nonlinear theory, virtually 
all the characteristics of the experimental pressure distribution are predicted qualita- 
tively by the theory (fig. 3). However, the effect of the entrainment is overpredicted 
somewhat. This error may arise from the failure to account for the curvature of the jet 
in computing the entrainment coefficients or it may be associated with an error in the 
empirical parameter a. In any case, if the entrainment coefficient is arbitrarily 
reduced by 40 percent, the correlation with experiment is much better for this particular 
example (fig. 4). The small discrepancy near the trailing edge on the lower surface 
results from the fact that the jet slot was located in the vicinity of the 98-percent station 
rather than at the trailing edge as assumed in the calculation. 

CONCLUDING REMARKS 

A nonlinear procedure for computing the pressure distribution on an airfoil with a 
trailing-edge jet flap has been described. The method is not restricted to thin airfoils 
or shallow jet-deflection angles. Correlation with experiment indicates that the 
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characteristics of the pressure distribution are predicted by the theory, but the effect of 
entrainment is overpredicted with the entrainment coefficient used. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 19, 1976 
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Figure 2.- Theoretical example demonstrating effects of supercirculation 
and entrainment on pressure distribution, r = 45°; Cj = 2.0; 
thickness ratio = 0.13. 
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Figure 4.- Correlation of experiment with theoretical pressures (modified 
entrainment coefficient). Thickness ratio = 0.125; Cj = 1.5; t = 31.4°. 
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